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Abstract 

We derive the exchange currents of pseudoscalar, vector, and scalar mesons from 
Feynman diagrams, and use them to calculate the magnetic form factors of nucleon 
and A(1232). The magnetic moments and electromagnetic radii are obtained by using 
those form factors and the parameters determined from the masses of nucleon and 
A(1232). We find the magnetic moments and electromagnetic radii of nucleon and 
A (1232) can be produced very well in the extended Goldstone-boson-exchange model 
(GBE) in which all of pseudoscalar, vector and scalar meson nonet are included. The 
magnetic moments of A(1232) are closer to experiment values and results from lattice 
calculation than the results obtained by the model without other mesons except for 
pion and sigma. 



1 Introduction 

In the frame of constitute quark model, the exchange currents between quarks are impor- 
tant to produce the properties of hadrons. There are two kinds of exchanges, gluon exchange 
and meson exchange. In Isgur's model, one-gluon-exchange (OGE) governs the structure of 
the hadronsp. Since the last few years, Shen et.aZ.[2J, Riska and Glozman [31 IH E] applied 
the quark-chiral coupling model to study the baryon structure. In the works of Glozman 
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et.al. [HE], the vector meson coupling was included to replace one-gluon-exchange. They 
found the spin-flavor interaction is important in explaining the energy of Roper resonance 
and got a comparatively good fit to the baryon spectra. Recently, L. R. Dai et. al. extended 
original chiral SU(3) quark model, in which the nonet pseudoscalar meson exchanges and 
nonet scalar meson exchanges are considered in describing medium and long range part of 
the interaction and OGE contribute the short repulse, to the extended chiral SU(3) quark 
model, in which vector mesons are included^. They calculated N-N scattering processes 
and found the similar results to the original chiral SU(3) quark model calculation. In the 1 Sq 
case, the one channel phase shifts of the extended chiral SU(3) quark model are obviously im- 
proved. In our previous works j7j , we found that mixing angles of N* resonance are sensitive 
to the different interaction models, and we calculated the amplitudes for photoproduction 
of negative parity N* resonances under 2GeV to compare different interaction models. All 
those suggest that vector meson exchanges may be important to the short range mechanism 
of quark-quark interaction. 

On the other hand, the method, by which the electromagnetic form factors of deuteron 
are calculated on nucleon level jH], has been applied to quark level by Buchmann et. aZ.[9 . 
They have calculated the electromagnetic property of nucleon and A(1232)|10 [ lllj. Their 
results of magnetic moments and electromagnetic radii of nucleon agree with experiment 
values well while the magnetic moment of A(1232) is 6.981/xjv- There are two experiment 
results, 6.14 ± 0.5lp N by LOPEZCATRO and 4.52 ± 0.50 ± 0A5p N by BOSSHARD [TJ. 
Through the result of Buchmann et.al. is still within the region given by Particle Data 
Group ^2], it is larger than both experiment values above. The recent result from lattice 
QCD is 4.99±0.56/iAr, which supposes the latter experiment result [I??] . In the calculations of 
Buchmann et.al., gluon, pion, scalar meson and confinement exchanges are considered. But 
the effects of other pseudoscalar mesons and vector mesons are not included. On nucleon 
level, the discrepancy between the coupling constant of ttNN (« 14) and that of pNN 
(~ 3) is large. The effect of p meson is negligible since there is square coupling constant 
in form factors. However, by using the relations: g nqq = and g pqq = (7pAw[14j, we 

know that the coupling constant of irqq is only about one fifth of ttNN while the coupling 
constant of pqq equals to that of pNN. Hence on the quark level, the effect of p meson 
on the form factors may not be ignored. In this paper, we calculate magnetic moments 
and electromagnetic radii of nucleon and A(1232) in the extended GBE model to study the 
effects of the vector mesons. 

In the following section, we present the Lagrangian we used. In section 3, the potential 
and masses of nucleon and A(1232) in the extended GBE will be calculated to determine the 
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parameters. The magnetic moments and electromagnetic radii will be obtained in section 4. 
Conclusion and discussion will be given in the last section. 

2 Lagrangian used 

First, we give the Lagrangians which will be used to derive potentials and the exchange 
currents. The total Lagrangian can be written as following: 

Ctotai = l -($Yd^ - d^Y^P) ~ m q W + E £<p + E £v + C*. (1) 

where ip and m q are the quark field and quark mass. Because we only calculate the case of 
nucleon and A(1232), we do not consider strange meson exchanges. 

Under a global infinitesimal chiral transformation, we can obtain the 7r part of Lagrangian 
as following: 

At = -g«$i E ^ • + ~ E D ^ a ■ DV > where : D ^ = d ^ a ~9pjl e abcP y.(2) 

a=l a=l b,c=l 

Here and henceforth Tc a and p° are the pion and rho meson fields, respectively, A" is the 
Gell-mann flavor matrix of i-th quark which reduce to Pauli matrix Tj as a=l,2,3 because 
we do not consider the strange dimension. g nq and g p are coupling constants. 

We construct the p quark coupling as Jido et. al. did in the case of p nucleon coupling 
[15J. According to the construction of the gauge theory, the direct coupling of a quark i[) 
to the p meson is constructed replacing the derivative by the associated covariant derivative 
with the gauge symmetry: 

d^(x) - D^{x) = d^(x) +ig P Y^{x\ V a {x)]pl, (3) 

a 

where V a is the generator of the gauge symmetry of the SU(3) flavor space. Here, [ip(x), V a (x)] = 
— ^\ a ip. Then the p part of Lagrangian is as following: 

C p = -g„tyf E A a • P a ^ - \ E P" ua ■ P%, where : p, v = d,p a v - d u p* -g.j^ e*4fi,(4) 

a=l 4 a=l b,c=l 

where we assume coupling constant between rho and quark, g pq = y , in order to be consistent 
to the 7r part in form. 

For rj and u mesons, Lagrangians have analogous forms as the part of ir and p, respec- 
tively. 
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We also consider a exchange : 

£ CT = 9a$o$ + ^(<9a) 2 . (5) 

The electromagnetic Lagrangians can be derived using the principle of minimal coupling 
dfj, — > d^+iqAfj,, where q is the charge carried by the field upon which the derivative operator 
acts. Then the 'jqq, '-/tttt and 7pp interaction Lagrangians are: 

£~fqq = -Qei^j^A^, (6) 
£ 77r7r = -ie[-Kd^ + - 7r + (9 M 7r]A M , (7) 



C ypp = -te[d^t(A^-A^p,)-d^p u (A^t-A,p^)] 



where 



7r = ^[n 1 +tn 2 }, p„ = \p\ + ip % ^ , (9) 

and Q is the charge of quark in unit of e. A^ is the electromagnetic field. Because the 
charges of ir°, p°, u, 4>, rj, r/and a are zero, there is no interaction between them and photon. 



3 The extended GBE model 

In the harmonic-oscillator model, we assume that the Hamiltonian of three quark system 
is of the form: 

# = E(™ 9 + ^-)+E vconf (^ ^ + E v res (^ (io) 

i=l Zm Q i<j i<j 

where are the spatial and momentum coordinates of the zth quark, respectively. The 

third term is a two-body harmonic-oscillator confinement potential: 

^(r l ,r,) = -a c ^-A J c (r l -r,) 2 , (11) 

and a>c presents the strength of the confinement and A^ is the color operator of the i-th quark. 
After removing the kinetic energy of the center of mass motion of three quark system, we 
take the first three terms of Eq. (JIUj) as unperturbed hamiltonian H . Hence the total baryon 
wave function &n(A) is an inner product of the orbital, spin-isospin, and color wave function 
and given by 

I ** ( A) >= (^T 3/2 exp(-(p 2 + A 2 )/26 2 )) | ST > N ^ x | [111] > N Jt\ (12) 
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where the Jacobi coordinates p and A are defined as p = (ri — r2)/\/2 and A = (rj + r2 — 
2r 3 )/v / 6. 

The residual interaction V™" 63 consists gloun and meson (pseudoscalar mesons tt, r], rj', 
vector mesons p, u, <p, and scalar meson a) exchange potentials. In the following, we will 
explain them in detail. 

The one-gluon exchange potential [16j: 



a,.„ ,J1 7r A 2 \ „, s 11 



4 2 J [ r m? q V 3 / 4m;? r 3 

1 1 



2m? r 3 



3 r x o(Pi - Pi) • o(°"i + °"i) - r x «(P* + Pi) ' ofo ~ 



,(13) 



where r = — r^; crj is the usual Pauli spin matrix. 
The one-7r exchange potential: 

V«(ij) = E K ■ Xfr ■ Vffj ■ V/ W (r), (14) 

0=1 

where 

/m v ) = . , \ 2 . (15) 

4TT{2m q ) z \ r r J 

M presents tt, 77, 77', p, u or <fi meson. /iM, and gMq is meson mass, cut-off and meson 
quark coupling constant. 

The one- 77(77') exchange potential can be obtained from one-7r exchange potential through 
replacing the all things about tt by corresponded ones of 77. 

The one-p exchange potential: 

V p °(ij) = E K ■ tyt x V • <r s x V/,(r), (16) 

0=1 

V p LS {ij) = E K ■ A, a 2[3r x \{p x - p 2 ) ■ i(crx + cr a ) - r x + p 2 ) • ~(<n - cr 2 )]/ p (r),(17) 

a=l ~ ~* 

V; C (u) = EAr-A^4m^(r). (18) 

a=l 

The one-o; and one-0 exchange potentials can be obtained from one-p exchange potential 
by analogous replacement. 

Here, we use the "perfect" mixing pattern for the 77, 77', u and <ft, i. e., 

\ 1 + V2, -1 + V2. x -1 + V2. 1 + V2 X nn , 
~ — ~^/6~ ~^ ^6 — ~~ V6 — ^ ~\/6~ 

A^ = W^A 8 + J^A , X( t> = ~yl Xs + yl Xo - (20) 



Hence we only consider nucleon, the strange dimension in the Gell-mann matrices above 
have no effect on the results. Then Eqs. (|19I20|) become: 



\ v — > — ^1, — T/^lj ^ — 1> ^<t> ~ 0? (21) 

that is, there is no contribution from </>. 

The scalar meson exchange potential is as following: 

a 2 ( ' p~^ r e~ Ar \ 

With the Hamiltonian of Eq. (|10j) and the wave function of Eq.([12|) it is straightforward 
to calculate the nucleon mass. One obtains 

M N {b) = Zm q + ^^ + V conf (b)-2a sS j^ + ^ 

M A (b) = ^ mq + 1 ^ + V conf {b)~2a s ^ + ^ 
where 

« 6) = (25) 

h ' (b) = -^{vtM^ 1 ~ )e" 2 -' 2 ' 2 erfc(^)) - ( m » Am)}, (26) 

where the individual terms in Eqs. (|23J24J) are the nonrelativistic kinetic energy, quadratic 
confinement, gluon, pion, eta, rho, omega and sigma contributions, respectively. In this 
paper we use rj to present both rj and r}'. Here 5 V presents + 8^). The confinement 
contribution to the nucleon and A mass is given by 

V conf (b) = 2Aa c b\ (27) 

and the a-meson potential contribution is 

v - {b) - -4vbs{( 1 - ^)^m^ 2 )) - « a)}. (28) 

In our calculation, we leave out all spin-orbit forces and the central components from the 
vector meson exchanges as Ref. jlj. This treatment is further supported on more theoretical 
grounds by a study of the two-pion exchange mechanism between constituent quarks in Ref. 

ini 
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Subtracting Eq. (j23j) from Eq. (j24|) all spin- independent terms drop out and one gets 

M A -M N = 6 g (b) + 5,(6) - l -8 v + 25 p - 6 U , (29) 

To obtain the numerical results, we show the parameters in our calculation. First we 
assume constituent quark mass m q = |Mjv = 313MeV. The gluon quark coupling constant 
a s is determined by Eq. (J2*9*|l . Here we use the same values of meson quark coupling constants 
for pseudoscalar and vector mesons respectively as Ref. j3], which is also used in Ref. jH]. 
The coupling constants between vector mesons and quark in Ref. j3] have included the part 
of tensor coupling. We follow them ont only in calculating mass but also in calculating 
magnetic moments and electromagnetic radii because the effect of tensor coupling part can 
be included by replacing vector coupling constants by the sum of vector and tensor coupling 
constant. This can be seen by comparing the exchange currents below of Gari jH] and ours. 
The values of Cut-offs A 7 follow the linear scaling lawsjl]: 

A 7 = A, + k (/i 7 — fi n ) for pesudoscalar mesons , (30) 
A 7 = A p + k (/i 7 — fi p ) for vector and scalar mesons . (31) 

The parameters b and a c are determined to get the right mass of nucleon. The explicit values 
of parameters are presented in Table ^ 

Table 1: Parameters: The quark mass is assumed to be one third of nucleon mass. The masses 
of mesons are from Particle Data group [T2]. The cut-offs, quark size and the coupling constants 



between mesons and quark 
nucleon. 


are from Ref. 


P]. b and 


a c are 


determined to get the right mass of 


m q 


313Mev 


Hit 


139Mev 


Hv 


547Mev /v 958Mev 


Hp 


770Mev 




782Mev 




680Mev 




0-67 [cf Vq 




1.31 




0.67 


K 


700MeV 




1200MeV 


K 


1.2 


a s 


0.8584 


b 


0.609fm 


a c 


17.98MeVfm- 2 



The contributions from each exchange potential is listed in Table El ot s here is smaller 
than Ref. ^Hj, 1.093. This is consistent with the conclusion of Ref. [Hj that the strength 
of OGE is reduced after the addition of vector mesons. This result a gap of contribution 
from gluon between this paper and Ref. (TOj about 120MeV. This gap is cancelled by the 
contribution from rho meson exchange. The contribution of cu-exchange is small. Then the 
mass of nucleon is reproduced with a small adjustments of b and a c . 
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Table 2: Contributions of the kinetic energy (without the rest mass term) and individual potential 
terms in the Hamiltonian to the nucleon mass of Ea. H23|) cc: color Coulomb part of V g , 5: ^-function 
part of V g . All entries are in [MeV]. 





r-pkin 


yconf 


V 9 

cc 


vi 


yTT 


yn 


V 




V a 


Total 


N 


503.1 


160.0 


-443.9 


39.6 


-127.2 


9.5 


-126.7 


24.8 


-39.4 


0.0 


A 


503.1 


160.0 


-443.9 


198.2 


-25.4 


-9.5 


-25.3 


-24.8 


-39.4 


293.0 



4 Magnetic form factors, moments and radii of nucleon 

By using the Lagrangians above, we can evaluate Feynman diagrams, which are shown 
by Figs. Al(a-e) in Appendix A, as Refs. [HI CHI HU EE] ■ The total current operator consists 
of the usual one-body operator and two-body exchange current operators tightly related to 
the different quark-quark interactions: 

3 3 

Ptotal^) = H Pimp(?i) + Y,{P 9 qq( r i' r j) + X! PMqq{*h *j) + Yl PSqq^i, Tj)}(32) 

i=l i<j M=ir,r],T]' ,p,u> S=cr,conf 

3 

J totalis) = 22 ^imp( r i) 

i=l 

3 

+ XX 3 gqq( r i' r i) + X JMqq(ri,rj)+ X ^MM^i, Tj) + £ 3 Sqg^i , Tj) .(33) 

i<j M=Tr,r],T]' ,p,u! M=n,p S=a,conf 

The explicit form of each current is given in Appendix A. 

Quite generally, the charge radius is defined as the slope of the charge form factor at 
zero-momentum transfer 

6 d 2 

i^(bW c(c 

where, according to the general definition of the elastic form factors jS], 



r c = -l^TTTo W) l q -o, (34) 



F c (q 2 ) = < JMj=J TM T | - 1 - / dft 5 p(q)y °(q) | JMj=J TM T > . (35) 

47T J 

We also consider the constituent quark size. It is an effective measure to including 
relativistic effects (TTJ [THj. We use (r 2 ) q = 0.36/m 2 [TT], which is consistent to Ref. [20J. 

Using Eq.(JH2J) and the ground state wave functions of Eq. (fT2"j) . we obtain charge radius 
of A + as following: 

r|+ =b 2 + r 2 g + -^(55 3 -S^-Sr,- 25 p - 25 w ) + -^V conf + r\. (36) 
n bm q brrig 



The charge radius of A is zero exactly, which is consistent with the analysis in Ref . JU] • 
If we compare this with the corresponding result for the proton 

4 = b2 + 4 + - S * ~ 2*,) + ^fconf + rl, (37) 

and neutron 

2 ^ (K j_ X 1 X x o; x N f2 MA-MN ,, fi , 

r n = - — (*, + ^ - -5, + 2* p - * w ) = -6 -^^r-, (38) 
we obtain from Eqs. (j36H3*S|) the parameter- independent result 

rl = r 2 v - r 2 n . (39) 

So we find that the relations Eqs. (38,39) found in Ref. JU] are still tenable through we 
consider both pseudoscalar and vector mesons. 

We list our numerical results in Table El The results agree with the experiment values 
well. By using Eq. (138(1 and experiment values, we can determined b=0.610fm, which agrees 
with the value b=0.609fm that we obtain in section 3 and use in the calculation. 



Table 3: Nucleon and A(1232) charge radii from individual two-body exchange currents. A finite 
electromagnetic quark size r 2 ^ = 0.36 fm 2 is used. The charge radius of the A is zero in the 
present model. The experiment values are from Particle Data Group|l2j. All entries are in [fm 2 ] 
except for the total result which is in [fm]. 





imp 


ggg 


~-.2 
nqq 


/v«2 
ijqq 


™-.2 
P<7<? 


™»2 
ujqq 


a 


/v«2 
conf 


\/\ r total 


Exp. 


V 


0.731 


0.094 


-0.060 


0.000 


-0.060 


0.000 


0.029 


-0.169 


0.751 


0.870 


n 


0.000 


-0.063 


-0.040 


0.007 


0.040 


0.020 


0.000 


0.000 


0.340 


0.341 


A 


0.731 


0.157 


-0.020 


-0.002 


-0.020 


-0.010 


0.029 


-0.169 


0.830 





The magnetic moments of A are defined as the q — > limit of the magnetic dipole form 
factor [S] 



= 2V67rMjV < JMj=J TM T | -1 / cM^q) x J(q)] 1 | JM J= J TM T >, (40) 
iq 47T J 



We calculate nucleon and A(1232) magnetic moments, and iV — > A transition magnetic 
moments. We consider next-to-leading order term of the isovector pion pair-current not only 
in the case of A(1232) but also in the case of nucleon, and find it is really no obvious effect on 
the result. The numerical results are presented in Tabled] It's found the magnetic moments 
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of nucleoli agree with experiment very well. The magnetic moments of A are proportional to 
their charge as only gluon, pion is considered. There are obvious improvements for magnetic 
moments of A. For A ++ , the effect of addition of 77 and vector mesons gives a decrease about 
1.3 from the values of Buchmann et. al., 6.9814// wjlUj. This value is more consistent with 
the experiment values, 6.14±0.51 by LOPEZCATRO and 4.52±0.50±0.45 by BOSSHARD 
[12J. It is just in the margin of the recent result from lattice QCD, 4.99 ± 0.56 (T^j. The 
values of magnetic moment for A + we obtained agree with the lattice values, 2.49 ± 0.47/xtv 
[12], and the experiment values, 2.7t\'^(stat.) ± 1.5(syst.) ± 3(theor) |2*T] compared with the 
value, 3.491, in Ref. |E|. 

There is only a small improvement for the N — > A transition magnetic moments. 
Our result is consistent with some authors, 2.44/ijv by simple nonrelativistic quark model 
(NRQM) [TJ] , 2.53(2.63)/^ by Dahiya et. aZ. |22J and 1.7 ~ 3.0/xjv by Julia-Diaz et. al. 
(2SI- The experiment values adopted by Dahiya et. al. is 3.1 /ii\r|2"4]. The recent one is 
3.642 ± 0.019 ± 0.085/xjv|2SI- The results of NRQM, Dahiya et. al, Julia-Diaz et. al. and 
ours are smaller than the experiment values especially the latter. The problem with the 
underestimation of the iV —>■ A transition magnetic moment persists also after inclusion of 
both pseudoscalar and vector meson exchange currents. 



Table 4: Nucleon and A(1232) magnetic moments, and N — ► A transition magnetic moments from 
individual two-body exchange currents. The experiment values of nucleon, A and N-A transition are 
fi p = 2.792847337(29)/i7v, n n = -1.91304272(45)/^, fi A ++ = 6.14 ± 0.51mjv (by LOPEZCATRO) 
and 4.52 ± 0.50 ± 0.45 (by BOSSHARD) pU, /t A + = 2.7t\i(stat.) ± 1.5(syst.) ± 3(theor) (21] and 
l L p^A+ = 3.642 ± 0.019 ± 0.085u^ |25j. respectively . All entries are in n^. 







f^gqq 


H"irqq 




t-^riqq 


f^pqq 




f^LLiqq 


M<7 




l^-total 


V 


3.000 


0.473 


-0.141 


0.469 


0.000 


0.000 


-0.051 


-0.148 


0.252 


-1.023 


2.830 


n 


-2.000 


-0.158 


0.195 


-0.469 


0.020 


-0.101 


0.051 


0.049 


-0.168 


0.682 


-1.897 


A++ 


6.000 


1.891 


0.325 


0.000 


0.121 


-0.604 


0.000 


-0.592 


0.504 


-2.045 


5.599 


A+ 


3.000 


0.945 


0.162 


0.000 


0.061 


-0.302 


0.000 


-0.296 


0.252 


-1.023 


2.799 


A 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


A" 


-3.000 


-0.945 


-0.162 


0.000 


-0.061 


0.302 


0.000 


0.296 


-0.252 


1.023 


-2.799 


p -> A+ 


2.828 


0.223 


-0.276 


0.663 


-0.057 


0.142 


-0.073 


-0.140 


0.237 


-0.964 


2.585 



The magnetic radius is defined as the slope of the magnetic form factor at zero momentum 
transfer: 

6 d 



r 2 



F M (0) rfq 2 
10 



q 2 =0, 



(41) 



we list our results of magnetic radii of the nucleon and A(1232) in Table |SJ The results 
agree with experiment very well. 

Table 5: Magnetic radii of the nucleon and A(1232) from individual two-body exchange currents. 
The magnetic radius of the A is zero. A finite electromagnetic quark size, r 2 q = 0.36 fm 2 , is used. 
Experiment values is from Refs. |26 | 127\ 12%]. All entries are in [fm 2 ], except for total results which 
are in [fm]. 





i gqq 


(>i2 
nqq 


(>i2 

77T7T 


r\qq 


/Vi2 

pqq 


■qpp 


2 2 2 

/y ^ ry ^ ty" 

uiqq a c 




rf 1 


Exp. 


p 


0.775 0.094 


-0.029 


0.217 


0.000 


0.000 


-0.011 - 


0.011 0.047 -0.331 





867 


0.855 


n 


0.771 0.047 


-0.263 


0.324 


-0.006 


0.011 


-0.016 - 


0.005 0.047 -0.329 





885 


0.873 


A 


0.783 0.191 


0.034 


0.000 


0.011 


0.220 


0.000 - 


0.022 0.048 -0.335 





830 






From Tables 3-5, we 


can find the effects 


of mesons 


except for pion and 


sig 


;ma 


is non 



negligible. They provide a contribution to give better result for magnetic moment of A. In 
other cases the result is similar to the one without those mesons since the effect of those 
mesons can be smeared by the variation of parameters. 

5 Conclusion and discussion 

We calculate magnetic moments and electromagnetic radii of nucleon and A(1232) in 
an extended GBE model. With the parameters determined by the mass of nucleon and 
A (1232), the magnetic moments and electromagnetic radii of nucleon are consistent with 
the experiment results very well. There are obvious improvements for magnetic moments 
of A. Our result is more consistent with experiment results and results from lattice QCD 
calculation compared with the model which only consider pion and sigma mesons. The 
relations obtained from one-pion exchange survive from addition of other mesons. The 
results above should be helpful in better understanding the short range mechanism of quark- 
quark interaction. 
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Appendix A: Charge density and exchange current 



n 



By using the Lagrangian in Eq. (1), we can evaluate Figs. Al(a-e) to obtain charge 
density and exchange currents as Refs. [TUl HU HHj.We obtain the charge density and 
exchange current as follows: 



(a) 



(b) 



(e) 



Figure Al: One-body and two-body exchange currents between quarks: (a) impulse, (b) gluon 
pair, (c) meson (tt, r/, r/', p, oj) pair, (d) mesonic, (e) scalar(confinement, a) pair. 

a) charge density 



Pimp ( r «> Q) 



Q, 



ee 



(Al) 



Pgqqi^ii ^ji Q) 



^ \1 ■ \ c 3 feee*™ [q • r + (<r, x q) • {* 3 x r)] + (i <- j. , r 



)}^( A2 ) 



Pirqqy^-ii ^jj Q.J ee 



-T ■ T J + T jz )Oi 



%q_ 

2m r 



fj ■ V r / W (r) + (z j), 



p w5 ( ri , r j; q) = ee iq <{- + r^cr; ■ • ^rf v {r) + (i j)> 



3' 



p pw -(r i; r i; q) = ee"^(^r J • r j + r^)^ x ^- ■ cr 3 x V P / p (r) + (z <-> j), 



P«5?ir», r,-, qj = ee 



/,1 



■'9 



2m„ 



o-j x V r / W (r) + (i <-> j). 



(A3) 
(A4) 
(A5) 
(A6) 



p 77r7r (r i; Tj, q) « 0, p 7PP (ri, r j; q) « 0. 



(A7) 



psO^r^q) = e^Q^e 



(2m,)3V5 q2 - iq ' V ' + 5 V 'l ,/s(r "^- (A8) 
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b) exchange current 



J imp (r,,q) = Qie-^-tifaxpue^ + ipite***}). (A9) 
2m q 

•W^r^q) = A< • X^ee^^a, + try) x r + (z ~ (A10) 







= ee i<J - ri ( 


A 


+ r i^) 4m 2 x Vr-cr i ' + ( ? 




(All) 


jirqqi^ii 




= ee* 9 '* 1 


[r l X T j 


]3<ri(<Tj • V r )K + (z ++j), 




(A12) 


jrjqqi^ii 




= ee iq - ri ( 




'4^2 x V^-V^+^j), 




(A13) 


jpqq( r i) 


r i,q) 


= ee iQXl 


A 


+ t iz )((t 1 + cr 2 ) x V r y p + (z <-> 


j), 


(A14) 


j pqq{*-ii 




= ee iqXl 


[t % X T j 


] 3 [V r - o-! x (cr 2 x V r ]/ P + (i * 




(A15) 


jwqq( r ii 


r i,q) 


= ee iqxi 


(~ + ^ 


)(cr 1 + cr 2 ) x V r / W + (z <-> j). 




(A16) 



2 L r —L/^r 

j^(Ti, r j7 q) = jjr^in x t 2 ] z <Ti ■ Vio-2 • V 2 J dve iq {R - vr \zJ-^- - z A ^j—)(A17) 



j 7PP (r i)ri ,q) = __^ [Tl xT 2 j 



q " M 

2 



[4m 2 , - o-! x Vi ■ o-2 x V 2 ] / due*^"^^ - z^-). (A18) 



Jf //y (r,,r i; q) = -^{Q^^, x q^^r,) + j)}. (A19) 

We have used the following abbreviations R = (r^ + r 3 -)/2, z m (q, r) = Lr + zwq, and 

«) = [±g 2 (l - 4t> 2 ) + ^u 2 ] 1 / 2 in Eqs. (A17, A18). 
In the above derivation, we use strong interaction form factors given by Adam et. al. 
j, and neglect the nonlocal terms as usual. 
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